Quantum particle displacement by a moving localized potential trap 
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We describe the dynamics of a bound state of an attractive 5-well under displacement of the 
potential. Exact analytical results are presented for the suddenly moved potential. Since this is 
a quantum system, only a fraction of the initially confined wavefunction remains confined to the 
moving potential. However, it is shown that besides the probability to remain confined to the moving 
barrier and the probability to remain in the initial position, there is also a certain probability for 
the particle to move at double speed. A quasi-classical interpretation for this effect is suggested. 
The temporal and spectral dynamics of each one of the scenarios is investigated. 

PACS numbers: 03.65.-w, 03.65.Nk, 03.65.Xp. 



I. INTRODUCTION 



Recent developments in nanotechnology allow displac- 
ing miniscule particles, which can be as small as an atom. 
These particles' relocation can be achieved either by op- 
tical tweezers [1| or by Scanning Tunneling Microscopy 
(STM) 01 . The STM moves an atom by creating a poten- 
tial well at its vicinity. The atom is then trapped in the 
tip of the STM's needle and can easily be relocated along 
with the tip's position (see Fig.l). Beautiful structures 
with incredible (sub angstrom) accuracy were achieved 
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Since the atom is a quantum particle, localization at 
finite space is always partial. The sudden activation of 
the trapping well could cause an atom loss like in an 
equivalent decay process 0, H, @. Moreover, in this pa- 
per we show that the sudden movement itself (not only 
the abrupt capturing) can be responsible for the atom es- 
cape. It is also shown that not only do some of the atoms 
remain (on the average) at their initial state, but some 
will move beyond the tip's influence at double velocity. 



Bound states subjected to sudden perturbations have 
been studied in relation to the so-called deuteron problem 
0. The tunnehng dynamics of a bound state has been 
reported in a time dependent well Q and after suddenly 
weakening the strength of the potential ^ . Here we de- 
scribe the transport of particles initially trapped in a well 
which is shifted at constant velocity along a waveguide. 
Under strong transverse confinement, the dynamics be- 
comes effectively one-dimensional whenever all relevant 
energies are much smaller than the excitation quantum 
in the radial direction. 
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FIG. 1: System schematic. A single atom displacement by an 
STM tip. 



II. THE MODEL 



To simplify the system, the well is modelled by a one- 
dimensional delta function potential well. It should be 
stressed that a ID negative(positive) delta potential is, 
in fact, an exponentially shallow potential well, and can 
model with great accuracy any well (barrier) whose phys- 
ical dimensions are smaller than the de-Broglie wave- 
length of the particle [1]. 

Prior to t = the well is localized at x = 0; however, 
for t > the well moves at constant velocity v. The 
potential well can then be formalized 
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Vix,t) 



-7(5 (x) for t <0 
-7(5 (x - vt) for i > 



(1) 



Thus, the system dynamics are fuUy characterized by 
the Schrodinger equation: 



,9V d'^ip 



dt 



dx^ 



+ V{x,t)ij 



(2) 



where we adopted the units h = 2m — 1. 
The dynamics begin with the bound state of the t < 
potential, i.e., 




V'(a;,t = 0) = Vr/Se 



-y\x\/2 



(3) 



III. FREE EVOLUTION OF THE BOUND 
STATE 

If it hadn't been for the potential well, the particle's 
probability density would spread out freely. If it is as- 
sumed that for t > the well is absent, and the Hamilto- 
nian becomes purely kinetic, then for t > the dynamics 
is free, and it can be obtained using the superposition 
principle 

/oo 
dx'Koix, t\x, t' = 0)i^{x', t' = 0) (4) 
-oo 

with the free propagator 



Ko{x, t\x',t') 



4T:i{t - t') 



2 ■ (x-rc')^ 

e Mt-t') 



(5) 



The time evolution of the initial bound state can be then 
simply written in terms of the Moshinsky function as 

^-0(2;,^) = ^/JJ2[^I {x,-ij/2,2t) + M {-x,-i-f/2,2t)]. 

(6) 

where the Moshinsky function reads [1, [l^, [lH 



e 2t 1 
M(x, k, t) := ?«(— z), z = — 
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in terms of the Faddeyeva function w{z), which is de- 
fined as w{z) := erfc(— iz). On physical groimds 
it is clear that each of the M functions corresponds to 
a freely time-evolved cut-off plane-wave. Such solution 
entails the well-known diffraction in time phenomenon, 
which consists of a set of oscillations in the density pro- 
file 

However, the imaginary wavector —ij/2 
makes such transients evanescent [l3| . leading to a uni- 
form expansion. 



FIG. 2: Transients developed in the density profile of a bound 
state \4>{x,t)\'^, for a moving delta with v = 1.5 and 7 = 0.7. 
The dotted line follows the trajectory of the well x — vt, while 
the dashed line corresponds to 2vt. Units h = 2m = 1 are 
used in all figures. 



IV. UNIFORMLY MOVING WELL 

For t > the propagator should be extended to the 
V ^ case, for which the corresponding propagator can 
be obtained using Duru's method and as well as by 
means of the path integral perturbation series [TH . It can 
be conveniently written in terms of the free propagator 
and a perturbation term, represented by a Moshinsky 
function 

K'g'\x,t\x',t') = Ko{x, t\x\t') 



+ (7/2)e 

xM {\x - vt\ + \x' - wt'l, -1-^7/2, t) (8) 

The time evolution of the state in Eq. ^ can be also 
obtained in close-form, using the integral [16l | 



dx'e'''''' 
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Vo-ik 



M{\x'\ + \x'\,-iVo,t) 
[M{\xlk,t)-Mi\xl^tVo,t)]. (9) 



Taking t' — and using Eqs. ([3]), ^ and Q, one can 
readily find 



^ivx /2 



-.^ (7/2)^ 



7 — iv/2 

X [M{\x - vt\,-v/2 - ^7/2, 2t) - M{\x - wt|, ^7/2, 2t)] 
7 + iv/2 

X [M{\x - vt\,v/2 - 47/2, 2t) - M{\x - uti, 17/2, 21)] , 

(10) 
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as the sum of a free term plus a perturbation. 

The adiabatic Massey parameter [l^, which distin- 
guishes the distinct dynamical regimes is therefore 
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(11) 



so that for 9 <^ 1 the adiabatic dynamics is recovered 
while 9^1 corresponds to the infinitely fast displace- 
ment of the well (free evolution). 

Indeed, the eigenstate of a moving delta well is [ij] 



-t^ivx/2^-'r\x-vt\/2 



(12) 



which for t — 0, becomes e™^^''^ijj{x, 0). The fraction that 
remains bounded for t ^ oo is 



l(V'I"^(0)|^(0))P = 



7 



16 



(7^ 



.,2^2 



(4 + 



2^2 ■ 



(13) 



for t;7 <C 1, the exponential becomes e™^/"^ ~ 1 in the 
spatial range of the initial bound state, and the overlap 
becomes unity. 



V. THREE SCENARIOS 

When t oo the three domains that were discussed at 
the introduction (the particles that remain at the vicinity 
of a; = 0, the ones that are localized to the well at uniform 
velocity and the ones that propagate at double velocity) 
eventually appear. The wavefunction can be written as 
a supperposition of three terms: 
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FIG. 3: Spatial distribution of the probability density 
(solid line) and its approximation (dashed line). The three 
domains (v = 0,v, and 2v) are marked. The parameters are 
7 = 10, « = 40 and t = 0.05. 



particle's velocity behave like Av ~ 7 and therefore, as 
can be seen in Eq. [141 the spatial width of the two peaks 
and 2v gets wider approximately like Aa; ^ jt (unlike 
the width of the localized part, which remains ~ 7~"'^)- 
Therefore, the distinction between the three parts can 
appear only when u > 7 (or 9 > 1). Moreover, due to 
their initial width, the peaks shape appears only when 
t > (■i;7)"^ 

The probability density of the exact solution with a 
comparison to the approximation, which focuses on the 
three terms is illustrated in Fig. 3. 



tp ~ tpfree + tpwell + ^2v 



(14) 



where 



1pfreeix,t) 
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injt 1 + {x/"ft) 



2 ' 
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\x — vt\ 



1 + w/27 {x - 2vt) + ijt 



(15) 



The first term ip f^ee describes the free evolution of the 
initial state in the absence of the well, as can be appre- 
ciated from the structure of the propagator. The second 
contribution remains localized in the moving trap and 
follows its classical trajectory x = vt, while the last term 
is responsible for the appearance of a peak in the density 
profile at X — 2vt. This contribution results from the 
partial reflection from the attractive well of the initial 
state probability density located at a; > vt. 

Since in the initial state, the particle was localized in 
a region as small as Aa; ~ 7^^ the uncertainty in the 



VI. ASYMPTOTICS 

The transition from the initial stationary bound state 
to the final moving one involves the two natural frequen- 
cies of the system: The frequency (energy) of the initial 
state /i — 7^/87r and the kinetic energy of the mov- 
ing particle /2 = u^/Stt. The a; = and the x = 2vt 
peaks are affected only by the frequency /2, however, the 
X — vt one oscillates with three harmonics /i, /2 — /i 
and 2/2 — fi- In Fig. 4 the temporal dynamics of the 
three peaks is shown. The x — and the x = 2vt decay 
like ^ i^^/^, while the x = vt one converges to its final 
constant value 



\4'{x = vt, t 



7/2 



(1 + («/27) 



2^2 ■ 



(16) 



In Fig. 5 a numerical spectral distribution ^'(/) = 
FFT[\ip\'^] (FFT stands for the Fast Fourier Transform) 
of each one of the peaks is presented. The four different 
frequencies are clearly shown. 

The 'A' peak corresponds to the frequency fi — j'^/Sn. 
The 'D' and 'E' peaks correspond to the frequency /2 = 
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FIG. 4: The temporal dynamics of the three peaks (a; = 0, 
X = vt and x = 2vt). The system's parameters are: 7 = 10, 
V = 20. 



It is also instructive to investigate the system in a mov- 
ing frame of reference, in which the well is at rest (origi- 
nally, at the lab reference the well moves to the right). 

Att = the particle (at the well's reference frame) be- 
gins to move to the left with respect to the well. We can 
regard it as three different scenarios: particle at the right 
of the well (gray) at the well (black) at its left (white) 
At t = all three types begin to moves simultaneously 
to the left at velocity v (Fig.6A). 

The white is free - so it remains at velocity v to the 
left. The black is trapped - so its average velocity is zero. 
But the gray hits the barrier and turns back with velocity 
—V, i.e., the final scenario is shown at Fig.6B. 

When we return to the lab frame of reference (where 
the well moves to the right), we see that the white one 
didn't move, the black moved with the well at velocity v 
and the gray moved with velocity 2v (Fig.6C). 

Obviously, this semiclassical interpretation is possible 
only due to the partial localization, which is a manifes- 
tation of Quantum mechanics. 
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FIG. 5: The spectral distribution *(/) of the three peaks (of 
Fig.4). 



uVStt, the 'B' one stands for h = I./2-/1I = Iu^-t^I/Stt 
and finally 'C stands for /4 = /g /a = \2v^ - 7^1 /Stt. 
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VII. SEMI-CLASSICAL REALIZATION 

Clearly, the double velocity effect cannot be classical, 
since a localized state will remain localized in the classical 
world. However, the origin of this effect has partially 
a classical interpretation. Quantum mechanically, the 
particle in the initial state is not completely localized 
inside the well. In fact, when the well is very narrow 
most of the chance is to find the localized particle outside 
the well. 



FIG. 6: Semiclassical realization of the three domains. 



VIII. CONCLUSION AND DISCUSSION 

We have presented a ID quantum model for an atom 
displacement with an STM tip. The model consists of 

a delta function well, which model the STM's potential 
at the vicinity of its tip end, which moves uniformly. It 
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was shown that the probabihty to remain trapped in the 



moving tip is 
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Moreover, it was shown that besides the trapped par- 
ticles, and the particles that remain close to their initial 
state, there is also a third group of particles, which prop- 
agates at double velocity (2u) away from their initial po- 
sition. We show that the probability for each one of the 
three groups has a different temporal dynamics. 
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